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Abstract. In this work we consider two models of two dimensional discrete 
systems subjected to three different types of coupling and analyse systematically the 
performance of each in realising synchronised states. We find that linear coupling 
effectively introduce control of chaos along with synchronisation, while synchronised 
chaotic states are possible with an additive parametric coupling scheme both being 
equally relevent for specific applications. The basin leading to synchronisation in 
the initial value plane and the choice of parameter values for synchronisation in the 
parameter plane are isolated in each case. 
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1. Introduction 

The presence of chaos has been extensively demonstrated in many naturally occurring 
nonlinear systems. The capability of the chaotic state to amplify small perturbations 
improves their utility for reaching specific desired states with very high flexibility and 
low energy cost. However there are many situations like power systems where control 
of chaos is desirable. This helps in improving a desired behaviour by making only small 
time dependant perturbations in a system parameter. In this context, a key observation 
is that a chaotic attractor has embedded within it an infinite set of unstable periodic 
orbits and the control mechanism aims at stabilising the system to any chosen one 
among them. The theoretical and experimental works so far done on this topic shows 
that various methods of control are capable of leading to improved performance and 
ordered response pQ. Moreover the naturally occurring chaotic orbits on a dynamical 
attractor has been utilised to carry information 0. Thus one can encode and transmit 
a prescribed message using chaotic states. This process may offer practical advantages 
over the usual periodic carriers, such as the possibility of real time reconstructions of 
signal dropouts in the communication. In such cases, it is desirable to have a number of 
systems behaving in a synchronised manner and this is achieved by coupling or feedback 
between systems. 

The phenomenon of synchronisation in chaotic systems has received much attention 
since 1990 [3 . It is found to occur in the dynamics of many coupled physical 
jU El El 13 El IE| and biological systems [101 HU H2] such as neurons in a network 
[T3*l IT^l 1151 ITS] . The synchronisation schemes form the basis of biological clocks that 
regulate daily and seasonal rhythms of living systems from bacteria to humans. Even 
when the individual systems are in the realm of chaotic behaviour, the coupled one 
most often can be synchronised in periodic or regular states. Thus synchronisation 
techniques automatically brings in control of chaos also. However there are applications 
like information processing where synchronised chaotic states are desirable. 

There are different types of synchronised states [13 HH] occuring in coupled chaotic 
systems. The first type called complete synchronisation has received much attention 
with the development of a new type of communication technique that exploits the 
possibility of masking a message by mixing it with a chaotic signal [Ej. With both the 
amplitudes and phases of the coupled systems varying in the same way the synchronised 
chaotic state will be restricted to a smooth invariant manifold of lower dimension 
than the full phase space. Complete synchronisation implies that the details of the two 
systems are identical or their difference converges to zero, pH E01 El E2j- However 
there are situations where one observes lag synchronisation in which the system states 
coincide only when one of them is shifted in time. So also there are cases of generalised 
partial synchronisation where no amplitude synchronisation is observed but the average 
amplitude is bounded and varies periodically [23J. 

The majority of the research works reported in the literature have been focused 
on synchronisation of continuous systems and a few on discrete systems concentrate on 
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one dimensional map |19| HU\ I2*T] . In this paper we consider two dimensional discrete 
systems and introduce different coupling schemes and the resulting possible types 
of synchronisation. Such systems model many interesting situations like population 
dynamics of mutually dependent species, neural networks, image processing techniques 
etc. Hence the possibility of their evolution in a synchronised manner under different 
types of connectivity forms a useful study. Moreover, two dimensional discrete systems 
occur as the Poincare maps of three dimensional continuous systems, which are the 
lowest dimensional systems capable of exhibiting chaos. Hence the synchronised 
behaviour of the latter can be analysed in a general and systematic way using the 
former. 

The paper is organized as follows. In § EJ we describe the different coupling schemes 
and discuss the stability of synchronisation for each one. In §01 and §|H these coupling 
schemes are applied to two specific examples of discrete systems viz. Gumowski-Mira 
map (GM) and Modified Gumowski-Mira map (MGM) maps and their synchronisation 
is studied using similarity function, stability analysis, synchronisation and stabilisation 
times and available basin in the parameter as well as initial value plane. It is found that 
complete synchronisation with control, lag synchronisation and synchronised chaotic 
states can be realised by choosing different coupling techniques. Concluding remarks 
are given in § 

2. Coupling schemes for 2-d maps 

We start with a two dimensional discrete system which is written as 

XI = F(Xl,n) andX2 = F(X2,fi). (1) 

where XI = (XI, Yl) define two dimensional space of the first map, X2 = (X2,Y2) 
the two dimensional space of the second map and /i is the control parameter. When 
both the systems are coupled the dynamics of the system develops as 

XI = F(X1, fj) + h (e, XT, X2) 

X2 = F(X2 lf i) + I 2 (e,XT,X2) (2) 

Here I\ or I 2 is the coupling scheme applied to both the maps which in general is a 
function of the variables of both maps and the coupling parameter e. Different types of 
/ can be chosen to achieve the required synchronised state. 

The possible forms for the coupling scheme / that we consider in this work are 

(i) linear coupling (CS1) applied to one of the equations, ie. 

h = eX j 

(ii) linear difference coupling (CS2), F = e(Xj — Xi) 

(iii) additive parametric coupling (CS3), F = fj,(l — eXj). 
For all the above cases i,j = 1,2 and i ^ j. 

The first two are well studied schemes of synchronisation while the third one is 
introduced for the particular type of two dimensional system discussed here. It is 
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interesting to note that this coupling depends on the parameter /i of the individual 
system in addition to the usual coupling parameter e. We find that this introduces 
more flexibility in the scheme by providing two parameters, one of which controls 
the dynamics also. This can result in a wider range in the (/i, e) plane for optimum 
synchronisation. The stability analysis is carried out and the similarity function plotted 
against time to characterise the nature and robustness of the synchronised state. In 
each case the efficiency of the scheme is quantified by computing the average time T\ 
taken for achieving synchronisation over a number of initial conditions. The robustness 
is checked by perturbing the synchronised state with a random noise and computing the 
stabilisation time T2 required to regain synchronisation after perturbation. 

Regions in the parameter plane (fi, e) that leads to different synchronised states 
are isolated and the basin of synchronisation in the phase plane marked out. These 
characterisations help to estimate the relative merits of different schemes. 

The stability of the symmetric state Xj = Xj and Y,i = Yj with respect to the 
transverse perturbations (5X, 5Y) is analysed by linearisation. We get 

SX1 = ^SXl + ^5Y1 + eC(5X2, 5X1) 



where C is in general a function of SX1, 5X2. Defining V = (5Xi — 5X 2 ) and Fij = J^. 



(i) For linear coupling this is obtained as 

V: +1 = J2 FijVi - eQ with C x = 1, C 2 = 0, % = 1, 2. (4) 

j 

(ii) For linear difference coupling, 

V: +1 = Y,F ll V^-2eC l (5) 

j 

with C\ = 1 and C2 = 0. 

(iii) For additive parametric coupling, 

K+l = Yl F ^ V n + ftCi With C t = 1,C*= I- (6) 
3 

Eigen values a 1 of the Jacobian matrix F^ are evaluated and the transverse 
Lyapunov exponents (AjJ calculated. 
For CS1, 

Al = lnK-e). (7) 

For CS2, 

X L = \n{a l -2e). (8) 

For CS3, 

Ai = huy + /ie). (9) 

The largest transverse Lyapunov exponent being less than zero indicates the stability 
of the synchronised state [H?] . 
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Figure 1. GM patterns generated in the phase space ie. (X,Y) plane by the 
discrete system in equation |[TDJ|. Here a = 0.008, b = 0.05, X = 0.1, Y = 0.0. 
—0.6 < fi < —0.1, A/i = 0.1. It is clear that the asymptotic dynamic states depends 
very sensitively on the control parameter \i. 



3. Synchronisation in Gumowski-Mira maps 

As a first example of a two dimensional discrete map we consider Gumowski-Mira Map 
(GM)|24J. This is a two dimensional difference equation in (X, Y) involving a non linear 
function f(X) that is advanced in time by one iteration in the Y equation. The iterates 
of this map give rise to a variety of interesting patterns in the (X, Y) plane and hence 
useful as a two dimensional iterative scheme for producing fractal objects. The GM 
transformation is defined in the (X, Y) plane as 

X n+1 =Y n + a(l - bY*)Y n + f(X n ) 

Y n+1 = - X n + f(X n+1 ) (10) 

with f(X n ) = »X n + 2{1 ~^f n . 

' n 

This map is found to have a very sensitive dependence on the control parameter 
fi I2H1I2E1- The phase space plots known as GM patterns in Fig. ((TK) justify this. In 
these calculations a and b are fixed to be 0.008 and 0.05 respectively and /i is increased 
in steps of 0.1. In an earlier analysis we have found that many recurring periodic and 
self similar sub structures, each with its own intermittency, periodicity, quasi periodic 
bands, band merging etc. in the bifurcation scenario of the system explain the 
dynamics of these patterns. 

Here we apply different coupling schemes mentioned in the previous section to the 
above system and try to target it to lower stable periodicities in their synchronised 
states. 
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Figure 2. Linear coupling applied to two GM maps a) The synchronised four cycle for 
H = —0.39. b) The points in the (//, e) plane that lead to synchronised four cycle, c) 
The error function e(X) ploted against iteration time n. When a random perturbation 
is applied after 6000 iterations synchronisation is regained in 410 steps. 

In the case of CS1 the pair of synchronising systems evolve as 
Xl n+1 = Yl n + o(l - bYl 2 n )Yl n + f(Xl n ) + eX2 n 

Yl n+1 = -Xl n + f(Xl n+1 ) (11) 
X2 n+1 = Y2 n + a(l - bY2 2 n )Y2 n + f{X2 n ) + eXl n 
Y2n+i — — X2 n + f(X2 n+ i) 

where e is the coupling parameter. The advantage of the method is that the feedback 
in the X equation alone suffices to achieve total synchronisation of both the variables. 
We start with \i = —0.39 (the individual systems are in the chaotic region), for e = 0.7 
and initial conditions (0.1, 0.0), (0.2, 0.12), the X and Y equations of both settle to 
an identical four cycle which is clear from Fig.(J2K)- Fig (J2t>) gives the (//, e) plot for 
synchronising values. The error function defined as e(X) = XI— X2 and e(Y) = Y\—Y2 
tends to zero when synchronisation is achieved j2ZI • We test this numerically by iterating 
equation (jllj) for values mentioned above and plotting e(X) against iteration number 
n (Fig (J2fc))- The synchronisation time T\ and stabilisation time Ti which is the time 
taken to stabilise after applying a random noise for 10 iterations are shown in the figure. 
The values of T\ and T2 are found for 10 different initial conditions and their average 
calculated. We get T\ = 778 and r 2 = 410. The largest transverse Lyapnunov exponent 
for the four cycle works out to be —0.1564. This indicates that this synchronised state 
is stable. 

However this coupling is active even after achieving synchronisation. Hence to make 
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Figure 3. Lag synchronisation observed in GM maps under Linear difference coupling 

a) (XI, X2) plot with iteration number n showing synchronisation with a lag of r = 4. 

b) Variation of similarity function with iteration number n for r = 4. c) Regions in the 
(/!, e) plane that can lead to lag synchronisation, d) Variation of similarity function 
as a function of coupling parameter e plotted after 10000 steps for fi = —0.23. Lag 
synchronisation is realised only for a range of 0.471 < e < 0.495. 



the scheme cost effective,we try a linear difference coupling (CS2) of the two systems, 

Xl n+1 = Yl n + o(l - bYl 2 n )Yl n + f(Xl n ) + e(X2 n - Xl n ) 

Yl n+l = - Xl n + f(Xl n+1 ) (12) 

X2 n+1 = Y2 n + a(l - bY2 2 n )Y2 n + f(X2 n ) + e(Xl n - X2 n ) 

Y2 n+ i = — X2 n + f(X2 n+ i) 

With fi = —0.23 both the systems are individually chaotic. When e = 0.48 and 
initial conditions (0.1, 0.0), (0.2, 0.12) they stabilise to an eight cycle, with a time lag 
r = 4 between them. Synchronisation achieved here is lag synchronisation which can 
be confirmed using similarity function, defined as (21 



Sf(r) = ([X2 " +T ~ X1 f 2 } (13) 
[(Xll)(X2l)f 2 

For a non zero value of r, Si(r) — > corresponds to lag synchronisation. It is clear from 
Fig (JHt)) the time series plot of (Xl n ,X2 n ) with r = 4. Fig (J2b) gives the Si(t) vs n 
plot which establishes lag synchronisation phenomena in the system. Fig (JSfc) gives the 
(fi, e) plot for synchronising values. Fig (jSJi) gives the variation of similarity function 
Si(t) with e. It is seen that lag synchronisation occurs only in a narrow region of e 
and the states are de-synchronised in the remaining regions. The largest transverse 
Lyapunov exponent, Aj_ = —0.431144, for the state with lag synchronisation. 
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Figure 4. Additive parametric coupling applied in two GM maps for \i = 0.001 a) 
e vs XI plot; the two system settles to four cycle for e = 0.8 indicating control of 
chaos before the onset of synchronisation; b) The synchronised periodic four cycle as 
e is increased to 0.85; c) The values in the (/i, e) plane that leads to synchronisation in 
the periodic four cycle; d) The evolution of the error function e{X) showing complete 
synchronisation and return to synchronisation after applying the random perturbation 
at the 7000 th step. 



The synchronisation in this scheme is only lag type. To achieve total 
synchronisation in a periodic state we try an additive parametric coupling (CS3). 

The dynamics of the coupled GM maps on applying CS3 develops through the 
equations 

Xl n+1 = Yl n + a(l - bYl 2 n )Yl n + f(Xl n ) + //(l - eX2 n ) 

Yl n+1 = - Xl n + f(Xl n+1 ) + /i(l - eY2 n ) (14) 

X2 n+1 = Y2 n + a(l - bY2 2 n )Y2 n + f(X2 n ) + //(l - eXl n ) 

Y2 n+l = -X2 n + f{X2 n+1 ) + fi(l - eYl n ). 

For /i = 0.001, both the maps individually exhibit irregular and aperiodic behaviour. 
After coupling with e = 0.85 and initial values (0.1,0.0), (0.2,0.12) they are found to 
settle to the same four cycle. Fig (jUi) gives the (e, x) plot for the first map which shows 
that for e = 0.8 the map settles to a four cycle ie. control of chaos is taking place before 
synchronisation and synchronisation is reached only when e is increased to 0.85. Fig 
(JUd) gives the (XI, X2) plot showing synchronisation and in fig (jit) the values of [i and 
e for which synchronisation is achieved are shown. 

The variation of the error function is calculated numerically by iterating equation 
(114}) and is plotted against iteration number n in Fig (@Ji). The average time taken for 
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Figure 5. The lag synchronisation observed under additive parametric coupling for 
e = 0.1, fi = —0.3 with an asymptotic 10 cycle (a) XI and X2 ploted with iteration 
number n showing a lag of r = 5; (b) the variation of the similarity function to 
zero after n = 1800 (c) the variation of similarity function as a function of coupling 
parameter e. Lag synchronisation is realised in the range 0.075 < e < 0.15. 



synchronisation is found to be T\ = 6265 and the average stabilisation time r 2 = 3844. 
The largest transverse Lyapnunov exponent is calculated as explained in the previous 
case and is obtained as Aj_ = —2.75098. Thus the synchronised state is a periodic four 
cycle which is stable. 

With n = —0.3, initial conditions (0.1,0), (0.2,0.12) and e = 0.1 both the systems 
are found to synchronise to a stable ten cycle. We note that this is also a case of lag 
synchronisation ie., Xl n = X2 n+T ,with r = 5. The time series plot of Xl n and X2 n in 
Fig (JSk) and Si(r) vs n plot in Fig (J^Jd) establishes this phenomenon in the system. We 
observe lag synchronisation for the range 0.075 < e < 0.15 (Fig (jEfc)). 

The range over which lag synchronisation exists in this case is larger than the 
previous type of coupling. However beyond e = 0.15, the two systems are totally 
desynchronised. 



4. Synchronisation in Modified Gumowski-Mira maps 

As a second example we consider a variant of the GM map which is hereafter called 
modified GM map (MGM map), defined by the following set of equations. 

X n+1 =Y n + a(l - bYl)Y n + f(X n ) 

Y n+1 = -X n + f(X n ) (15) 
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Figure 6. Linear coupling applied to two MGM maps a) Regions in (/j,, e) plane leading 
to synchronised states, b) basin of the synchronised states in the initial value plane, 
c) error function e(X) vs iteration number n plot for e = 0.4 where the synchronised 
state is a fixed point. A random noise is applied at the 6000 th iterate. Synchronisation 
is regained after 98 iterations, d) error function e(X) vs n plot for e = 0.5 where the 
synchronised state is a limit cycle. A random noise is applied at the 6000 th iterate. 
Synchronisation is regained after 520 iterations. 



with f\X n )=^X n + 2{l ^f" . 

' n 

In this map the iterations in both X and Y are in step. With values of a and b are fixed 
as a = 0.008, b = 0.05, the system has a prominent period doubling sequence in the 
bifurcation scenario In the GM map considered in § El above synchronised states 
are periodic states which means synchronisation follows control of chaos. We find that 
synchronised chaotic states can be realised in addition to synchronised fixed points and 
limit cycles in the modified GM map introduced here. 

We apply CS1 to two such systems. In this case for /i = —0.7 and initial conditions 
(0.1, 0), (0.2, 0.1) both the maps are individually chaotic. With coupling for e = 0.4 they 
asymptotically approach a fixed point (1.66, —0.3316), while for e = 0.5 they synchronise 
to a limit cycle. 

Fig (JB^) gives the (/U, e) plot where synchronisation sets in. The region to the left 
of the white line corresponds to synchronisation to a fixed point while that to the right 
corresponds to limit cycle. Fig (jHb) gives the basin of synchronisation of two MGM 
maps for e = 0.5 (when the final synchronised state is a limit cycle). Fig (jHt) and 
Fig (jntO gives the e(X) — n plot of coupled MGM maps from where synchronisation 
time and stabilisation time can be read off. For fx = —0.7 and e = 0.4 the largest 
transverse Lyapunov exponent (AjJ is (—1.883383) where one cycle periodicity is seen 
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Figure 7. Linear difference coupling applied to two MGM maps, a) (XI, X2) plot 
showing synchronised chaotic states b) regions in the (fx, e) plane where synchronisation 
occurs, c) Variation of the error function e(X) with n for li — —0.7. The time taken to 
reach synchronisation and the time taken to regain the same after a noisy perturbation 
are more in this case d) basin of synchronised chaotic state in the initial value plane. 



in the synchronised state. For e = 0.5 when the synchronised state is a limit cycle the 
value A ± is (-1.44478). 

We try CS2 for two MGM maps. For fi = —0.7 both the systems are individually 
chaotic and for e = 0.85 we observe chaotic states synchronised in the coupled system. 
This is shown in (XI, X2) plot in Fig (jZk). Fig (J2b) gives the regions in the (//, e) plane 
leading to this type of synchronisation. Error function e(X) and e(Y) evolve into a 
stable fixed point (0,0), which is shown in Fig fZt). The average synchronisation time 
and stabilisation time for 10 different initial conditions are found to be T\ = 2605 
and Ti = 1057. Fig gives the basin for synchronisation in the range [-5,5]. 

Aj_ = —0.40786 shows that the chaotic synchronised state is stable. 

We try CS3 scheme to two MGM maps. 

Here for \i = —0.3 initial conditions (0.1, 0.0)(0.2, 0.1) (where individual systems 
show 4 cycle periodicity) for e = 0.5, we observe synchronised chaotic states in the 
coupled systems. This is shown in (XI, X2) plot in Fig (jH^i). Fig (jHjo) gives the regions 
of the parameter plane (fi, e) leading to this type of synchronisation. 

e(X) evolve into a stable fixed point (0,0) which is shown in Fig (JBfc). A similar 
plot is obtained for e(Y) also. Time taken for synchronisation t\ = 282. After 2000 
iterations a noise is applied and stabilisation time is obtained as T2 = 186. Fig (JEtl) 
gives the basin for synchronisation in the range [-5,5]. 

For this coupling, the largest transverse Lyapunov exponent computed works out 
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Figure 8. Additive parametric coupling in two MGM maps a) The (XI, X2) plot of 
the synchronised chaotic state b) points in the (fi, e) plane resulting in synchronisation, 
c) Variation of error function e(X) with time for /i = —0.3. d) The basin of 
synchronised chaotic state in the initial value plane. This scheme is more efficient 
and the state is robust to perturbation. 



to be —1.227154 indicating stability for synchronised chaotic state. 



5. Conclusion 



In this work, a study on different coupling schemes for synchronisation in two 
dimensional systems is carried out. We consider two types of two dimensional 
discrete maps called Gumowski-Mira(GM) maps and modified Gumowski-Mira(MGM) 
maps. Three different coupling schemes namely linear coupling(CSl), linear difference 
coupling(CS2) and additive parametric coupling(CS3) are applied to them. In GM maps 
using CS1 and CS3 control of chaos is achieved along with complete synchronisation. 
However only lag synchronisation is observed with CS2. In MGM maps while CS1 
stabilises the maps to fixed points or limit cycles, CS2 and CS3 lead to synchronised 
chaotic states. It is found that the synchronisation time and stabilisation time is shortest 
for CS1 in GM maps compared to CS3. In the case of MGM maps also CS1 offers the 
shortest synchronisation time and stabilisation time. In this case both CS2 and CS3 
yields chaotic synchronised states among which T\ and T2 are shorter for CS3. It is found 
that the basin of synchronisation in the initial value plane is unique for all the three 
types of coupling in GM maps, while the regions leading to synchronisation in MGM 
maps are small but finite. The parameter values (n, e) leading to synchronisation in GM 
maps has largest range for CS3 while those leading to lag synchronisation in CS2 has 
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limited range. In MGM maps this range is maximum for CS1 and minimum for CS2. It 
is found that although CS2 is cost effective, CS1 gives maximum range in (//, e) plane, 
wider basin in (AO, YO) plane and shorter times for synchronisation and stabilisation 
after perturbation. CS3 has wider range in parameter with control of chaos in GM map 
and is useful to get synchronised chaotic states in MGM maps. 

Although control of chaos is simultaneous with synchronisation in most cases, 
synchronised chaotic states are possible in the second example considered here. This 
may find applications in information processing using connected systems. 
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